The land of deformation south of 68 Ni 
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We study the development of collectivity in the neutron-rich nuclei around N = 40, where exper- 
imental and theoretical evidences suggest a rapid shape change from the spherical to the rotational 
regime, in analogy to what happens at the island of inversion surrounding 31 Na. Theoretical calcu- 
lations are performed within the interacting shell model framework in a large valence space, based 
on a 48 Ca core which encompasses the full pf shell for the protons and the O/5/2, IP3/2, lPi/21 Ogg/2 
and lds/2 orbits for the neutrons. The effective interaction is based on a G- matrix obtained from a 
realistic nucleon-nucleon potential whose monopole part is corrected empirically to produce effective 
single particle energies compatible with the experimental data. We find a good agreement between 
the theoretical results and the available experimental data. We predict the onset of deformation at 
different neutron numbers for the various isotopic chains. The maximum collectivity occurs in the 
chromium isotopes, where the large deformation regime starts already at N = 38. The shell evolu- 
tion responsible for the observed shape changes is discussed in detail, in parallel to the situation in 
the N = 20 region. 

PACS numbers: 21.60.Cs, 21.10.-k, 21.10.Rc 
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I. INTRODUCTION 

In the last decades, more and more experimental evi- 
dence has been accumulated establishing the breaking of 
the shell closures known at the stability valley when ap- 
proaching the drip lines, mainly at the very neutron-rich 
side. The first example of an unexpected disappearance 
of a shell closure (N = 8) was found in 11 Be, whose 
ground state is an intruder 1/2+ located 320 keV below 
the "natural" Ofiw 1/2" state Q. However, the true rel- 
evance of this finding was not recognized till many years 
later, and even now it is shadowed by the fame of its more 
neutron-rich isobar 11 Li. In the sd shell, the expected 
N = 20 semi-magic neutron rich nuclei turned out to be 
actually well deformed 0, Q . The presence of deformed 
ground states in the nuclei around 32 Mg and the break- 
ing of the N = 20 shell closure in this region has been 
extensively studied experimentally and theoretically (see 
for shell model reviews) . The aim of these studies 
is to map the limits of the so-called island of inversion, 
i.e. the region of nuclei where the strong quadrupolc 
correlations overcome the spherical mean-field gaps, fa- 
voring energetically the deformed intruders, which often 
become ground states. Indeed, another basic aim is to 
understand microscopically the dynamics responsible for 
these shape transitions. 

The question of the persistence of the N = 40 har- 
monic oscillator closure in neutron-rich nuclei comes in 
naturally in this context. The high lying 2 + state ob- 
served in 68 Ni and its low B(E2; 2+ — ► + ) value are 
the result of the relatively large energy gap separating 
the pf and 0g 9 / 2 orbitals @|. However, this gap gets 
reduced (or even disappears) when protons are removed 
from 68 Ni: The nucleus 66 Fe, with only two protons less, 



shows a sudden change in nuclear structure with an in- 
creased collectivity manifested via its very low lying 2 + 
state. Along the iron chain, indications for a collective 
behavior come from the systematics of the 2 + states Q 
as well as from the recent measurement of the B(E2) 
values in 64 < 66 Fe [1, d|. The evolution of the B(E2) 
values in iron isotopes points to a sudden increase of col- 
lectivity when approaching N — 40. Only very recently 
the first measurement of the excited levels in 64 Cr has 
been reported [ll|. This is the nucleus where theoret- 
ical calculations predicted the lowest lying 2 + level in 
the region [l2l - fl4| . The measured 2 + state energy agrees 
within 100 keV with those theoretical predictions. The 
sudden shape change at N = 40 challenges the theo- 
retical models, which have to account for a particularly 
rapid shell evolution responsible for these effects. Re- 
cent beyond mean-field HFB+GCM calculations with the 
Gogny force, reported in Ref. [lj|, show an increase of 
collectivity towards the proton drip line. The spherical 
neutron single-particle energies obtained in the HFB ap- 
proach reveal almost no variation of the N = 40 gap 
with the proton number between Z = 20 and Z = 28. 
As a consequence, only a moderate collectivity in the iron 
and chromium chains is found, without visible structure 
changes between them. The collectivity at N = 40 has 
been previously a subject of many shell model studies us- 
ing different valence spaces and interactions 
In particular, it has been shown that the shell model cal- 
culations using the Of 5/2, l?>3/2; I.P1/2, O59/2 neutron or- 
bits (the fpg valence space) and realistic interactions [l3| , 
can reproduce rather well the level schemes of 62,64 Fe, but 
fail to do so for the 2 + state of 66 Fe [16j . To reproduce 
the large quadrupole collectivity in this mass region, the 
inclusion of the neutron le/5/2 orbital is needed, as first 
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surmised in Ref. [l2[ and confirmed recently in Ref. 
This can be explained in terms of the quasi-SU3 approx- 
imate symmetry: In this framework, the deformation is 
generated by the interplay between the quadrupole force 
and the central field in the subspace consisting on the 
lowest Aj = 2 orbitals of a major shell (T3] ■ 

In this work we discuss in detail how the sudden on- 
set of collectivity can be interpreted in terms of shell 
model calculations in large model spaces and we look for 
the similarities in the deformation-driving mechanisms 
at N = 20 and N = 40 when approaching the neutron 
drip line. We present novel calculations for the iV = 40 
region, in a model space including the pj-shcll for pro- 
tons and the 0/ 5/2 , lp 3 / 2 , lpi/2, 0g 9/2 and ld 5/2 orbits 
for neutrons. We discuss first our model space and in- 
teraction as well as the computational details in Sec. HU 
The shape evolution along the TV = 40 line is studied in 
Section Hm Next, we illustrate the development of defor- 
mation along the isotopic chains of iron and chromium 
in Sec. IIVI and we discuss the structure of the nickel 
isotopes. Conclusions are given in Sec. |V] 



II. THE VALENCE SPACE 

The physics of the N = 40 nuclei has been discussed 
in the recent past both experimentally and theoretically 
d 0, H Effig. In 68 Ni, the N = 40 harmonic oscillator 
shell closure is somewhat weakened, however, whereas 
some of its properties seem consistent with a supcrfluid 
behavior, others may point to a double magic character. 
Below 68 Ni, the iron and chromium isotopic chains are 
characterized by an open proton shell which favors the 
development of quadrupole correlations. Such a situation 
occurs also for the N = 20 nuclei in the so called island of 
inversion around 32 Mg fl8j : the protons occupy the N = 
2 harmonic oscillator shell and the neutrons may lie cither 
in the N = 2 harmonic oscillator shell (normal filling) or 
in the lp3/ 2 and O/7/2 orbitals from the N = 3 major 
shell (intruder configurations). The latter orbits form 
a quasi-SU3 block |19f which enhances the quadrupole 
correlations. Extended shell model calculations by the 
Tokyo group, with the SDPF-M effective interaction fl8| . 
predict the dominance of different type of configurations 
(OpOh, lplh, 2p2h) for the ground states of the nuclei 
with 18 < TV < 22 and 10 < Z < 14, delineating the 
contours of the island of inversion. 

To assess if such a scenario could develop around 68 Ni, 
we adopt here a model space based on a 48 Ca core which 
comprises the p/-shcll for protons and the lp3/ 2 , lpi/ 2 , 
O/5/21 0<?9/2) and le?5/ 2 orbits for neutrons. The degrees 
of freedom that can be encompassed in this valence space 
are very similar to the ones present in the study recalled 
previously [18j . and the same kind of phenomena could 
therefore be described. The advantage of this valence 
space is that it contains all the physical degrees of free- 
dom important for the description of the low-lying prop- 
erties of these nuclei, it is computationally tractable, (see 



precisions hereafter) and it is almost free of center of mass 
spuriousity since its main components, the Of 7/2 0g$/ 2 
excitations, arc excluded from the space. 

The interaction proposed in this work, denoted here- 
after LNPS, is a hybrid one, based on several sets of real- 
istic two-body matrix elements (TBME) . Its main build- 
ing blocks are: 

• The last evolution of the Kuo-Brown interaction 
(KB3gr) for the p/-shell 0; 

• The renormalizcd G- matrix of Ref. 1 2 ill with the 
monopole corrections introduced in |22j , for the re- 
maining matrix elements involving the lp3/ 2 , lj?i/ 2 , 
O/5/2 j and 0<7g/ 2 neutron orbits. 

• The G-matrix based on the Kahana-Lee-Scott po- 
tential [23|, for the matrix elements involving the 
1^5/2 orbit. This potential has been successfully 
employed in the definition of the recent SDPF-U 
shell model interaction [24[ for the description of 
neutron-rich sd — pf nuclei. 

In addition, another set of experimental constraints 
has been taken into account for the final tuning of the 
monopole Hamiltonian: 

• The Z = 28 proton gap around 78 Ni is inferred 
from recent experimental data in 80 Zn and fixed so 
as to reproduce the measured B(E2). A standard 
polarization charge of 0.5e has been used in all the 
calculations presented in this work. 

• The size of the N = 50 neutron gap in 78 Ni has been 
estimated to be ~ 5 MeV. The evolution of the neu- 
tron gap with the neutron number is rather inde- 
pendent on the pro ton number. On the contrary, 
the systematics [25( shows that the Ogg/2 — 1^5/2 
gap in Zr isotopes increases by 3 MeV when the 
neutron Ogg/2 orbital is filled. We assume a similar 
behavior for the Ni chain. With this assumption, 
the observed 5 + , 6 + states in 82 Ge, which are sup- 
posed to be lp-lh excitations across the N = 50 
gap, are correctly reproduced [26| . 

Finally, in order to compensate for the absence in our 
space of the third member of the quasi-SU3 sequence, the 
2s\/2 neutron orbit, we have increased the quadrupolc- 
quadrupole interaction of the neutron 0*79/2 and 1^5/2 
orbits by 20%. 

The effective single particle energies (ESPE) of the 
LNPS interaction are shown in Fig. [T]for the neutron or- 
bits at N = 40, between 60 Ca and 72 Ge. For comparison, 
we show in the same figure the neutron ESPE at N = 20 
between 28 O and 36 S, calculated with the SDPF-U inter- 
action [24|. The similarities are striking. In the N = 20 
case, a reduction of the neutron 0g?3/ 2 — O/7/2 gap takes 
place when protons are removed from the proton 0c? 5 / 2 
orbital. This feature, accompanied by the proximity of 
the quadrupole partner neutron orbitals Of 7/2 and lj?3/ 2 , 
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FIG. 1: (Color online) Neutron effective single particle ener- 
gies obtained with the LNPS interaction at N = 40 (upper- 
part) and with the SDPF-U interaction at N = 20 (lower 
part). 



is responsible for the formation of the island of inversion 
at N — 20. At N = 40 one observes the same behavior 
but for the neutron O/5 12 — Ogg 12 gap when going down 
from 68 Ni and the closeness of the quadrupole partners 

099/2 ~ 1^5/2- 

The analogy in the shell evolution between N = 20 
and N = 40 suggests therefore a possibility of another 
island of inversion below 68 Ni. 

As mentioned already, the deformation driving role of 
the neutron ld 5 / 2 orbital below 68 Ni was already dis- 
cussed in the shell model calculations of Ref. |12j . How- 
ever, due to the computational limitations at that time, 
the calculations were performed with a closed neutron 
lps/2 orbit, therefore using 52 Ca as a core. This is no 
longer an issue here: The present calculations have been 
carried out including up to 14 particle-14 hole excitations 



across the Z — 28 and N = 40 gaps, when it appeared 
necessary to assure the convergence of the calculated elec- 
tromagnetic properties. The largest dimensions of the 
matrices treated here reach 10 10 in the case of 64 Fe. All 
the calculations of this work have been performed using 
the m-scheme shell model code ANTOINE (13. 



III. THE PROPERTIES OF THE N = 40 
ISOTONES FOR Z<28 AND THE NEW REGION 
OF DEFORMATION 

The change of structure in even-even N = 40 isotones 
from 68 Ni down to 60 Ca is illustrated in Fig. [2] In part 
(a) the evolution of the excitation energy of the 2 + state 
is plotted while part (b) shows the corresponding evolu- 
tion of the E2 transition rates. The available experimen- 
tal data are satisfactorily reproduced in all cases. Our 
theoretical approach predicts the maximal deformation 
at the middle of the proton shell, i.e. in 64 Cr, where the 
calculated 2 + energy is the lowest and the value of the 
B(E2;2+ -> 0+) is the largest. 

In Table [I] we list the extra occupancies of the two 
neutron intruder orbitals, Ogg/2 and ltis/2, relative to the 
normal filling. The occupation of the Ogg/2 orbit in 68 Ni 
is close to 1, which corresponds to «50% of the calculated 
wave function having a J = pair excited across N = 40 
gap. The wave function is however difficult to interpret, 
because, in spite of its very large content of excited pairs, 
the doubly magic component is still substantial. This is 
probably the reason why 68 Ni shows at the same time an 
increase of the 2 + excitation energy typical for doubly 
magic nuclei, and no sign of shell closure in the neutron 
separation energy. 

The occupation of both neutron intruder orbits grows 
rapidly when protons are removed, due to the reduction 
of the neutron N = 40 gap shown in Fig. [TJ The ratio 
of the Ogg/2 to ld§/2 occupations also evolves from Ni to 
Ca, partly as an effect of the level crossing taking place 
around 62 Ti, but mostly due to the quasi-SU3 structure 
of the intruder states. The very important role of the 
ld 5 /2 neutron orbit in the build up of collecti vity in this 
region sheds doubts about the results of Ref. [Tj| , which 
do not include this orbit. The more so when one examines 
the very unrealistic ESPE that they enforce into their 
schematic interaction in order to get results close to the 
experimental data. 

In Table [J we list as well the percentages of neutron 
n-particle n-hole excitations in the ground state wave 
functions. The 4p-4h components are dominant in Fe, 
Cr, Ti and Ca, however 2p-2h and 6p-6h contributions 
are sizable. The complexity of the wave functions consti- 
tutes the main difference between N = 40 and N = 20 
regions. In the latter, nearly pure 2p-2h components has 
been shown to dominate the + ground states [l8[ . 

It should be also pointed out that this evolution of 
the neutron filling and of the particle-hole structure does 
not mean that the nuclei will become more and more 
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IV. THE EVOLUTION OF THE DEFORMATION 
ALONG THE ISOTOPIC CHAINS 

The isotope chains provide a very illustrative picture 
of the evolution of deformation in this region. In Fig. 
[3] we show the results obtained for the iron chain. In 
the panel (a) the excitation energies of the 2 + states 
are compared to the available experimental data. In 
panel (b) we show the ratio of the excitation energies 
E(4 + )/E(2 + ) and E(6 + )/E(i + ), which make it possi- 
ble to recognize whether a nuclear spectrum is close to 
that of the rigid rotor or not. Let us remind that the 
perfect rotational limit would require i?(4 + )/£'(2 + )=3.33 
and E(6+)/E(A+)=2.1. Panel (c) depicts the theoretical 
and experimental B(E2) transition probabilities along 
the yrast bands. The measured B(E2; 2+ — > + ) val- 
ues from Ref. @ are plotted with dots and from Ref. [To| 
with crosses. Finally, the intrinsic quadrupolc moments 
derived from the calculated spectroscopic ones are shown 
in part (d). To establish a connection between the labo- 
ratory and the intrinsic frames we use the relations: 



(J+l)(2J + 3) 
3K 2 - J(J+ 1)' 



3 (J), K?l, (1) 



and 



FIG. 2: (Color online) Evolution of the nuclear structure 
along the N = 40 isotonic chain: (a) theoretical 2 + exci- 
tation energies and (b) B(E2; 2 + — > + ) values, compared to 
the available experimental data. 



TABLE I: Occupation of the neutron intruder orbitals and 
percentage of particle-hole excitations across the N = 40 gap 
in the ground states of the N = 40 isotones. The last column 
contains the correlation energies evaluated for these states. 

Nucleus vgg/2 vd 5 /2 OpOh 2p2h 4p4h 6p6h E corr 



b8 Ni 


0.98 


0.10 


55.5 


35.5 


8.5 


0.5 


-9.03 


66 Fe 


3.17 


0.46 


1 


19 


72 


8 


-23.96 


64 Cr 


3.41 


0.76 





9 


73 


18 


-24.83 


62 Ti 


3.17 


1.09 


1 


14 


63 


22 


-19.62 


60 Ca 


2.55 


1.52 


1 


18 


59 


22 


-12.09 



deformed with decreasing Z] the ground-state deforma- 
tion properties result from the total balance between 
the monopole and the correlation energies (mainly of a 
proton- neutron character). In Table HI we list also these 
correlation energies extracted from the multipolc Hamil- 
tonian along the N = 40 line. Indeed, the correlation 
energy, reflecting the deformation, increases from Ni to 
Cr, where it reaches its maximum, and then diminishes 
toward Ca. The transition between Ni and Cr is not 
gradual: The removal of two protons already provokes 
and abrupt change, from spherical to a strongly deformed 
prolate shape. 



B(E2, J -> J - 2) = ^e 2 \(JK20\J-2,K)\ 2 Q 2 mt (2) 
for i, 1. 

The same properties are plotted in Fig. [4] for the 
chromium chain. An excellent agreement with the ex- 
periment is found for all excitation energies. 

The known B(E2) transition rates are well reproduced 
within the error bars as well. Comparing the results for 
both isotopic chains it can be seen that the onset of de- 
formation occurs at different neutron number. In the 
chromium chain the intrinsic quadrupole moment stays 
constant along the yrast the band already at N = 38. 
This is one of the fingerprints of a good rotor behav- 
ior. The iron isotopes undergo the transition at N = 40. 
We have also verified that the intrinsic quadrupole mo- 
ments obtained from the spectroscopic ones (Eq. (JTJ) are 
nearly equal to those obtained from the transition prob- 
abilities according to Eq. © in 6 2 .64,66 Cr and 66,68 Fe _ 
We obtain a value of Q mt ~ 150efm 2 in 62,64,66 Cr which 
corresponds to /3 ~ 0.35. The deformation of 66 ' 68 Fe is 
slightly lower, with Qi n t = 145efm 2 (ft ~ 0.3). In the 
lighter Fe and Cr nuclei, the discrepancies between the 
calculated quadrupole moments and those obtained in 
the rotational scheme from the B(E2)'s are much larger. 

Let us further discuss the structure of the calculated 
states. In Table [TT| we show the occupancies of the 
neutron intruder orbitals in the ground states of the 
chromium and iron chains. In the case of the 0c/g/ 2 or- 
bit the extra occupancy is reported, i.e. the difference 
between the value obtained in the configuration mixing 
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FIG. 3: (Color online) Theoretical results along the iron iso- 
topic chain in comparison with the available experimental 
data: the excitation energies of the 2 + states are shown in 
panel (a), in panel (b) we present the ratio of energies of 
E(J+2)/E(J), the B(E2) transition rates are plotted in panel 

(c) and the calculated intrinsic quadrupole moments in panel 

(d) . Two experimental sets of the B(E2) values are shown: 
from Ref. 0] (black dots) and from Ref. (T^] (crosses). 



FIG. 4: (Color online) The same as in Fig. [3] but for the 
chromium isotopes. 
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TABLE II: Extra occupations of neutron Ogg/2 an d W5/2 or - 
bitals in the ground states of the chromium and iron chains. 



Nucleus TV ^ld 5 / 2 



b2 Fe 


36 


0.95 


0.12 


64 Fe 


38 


2.0 


0.27 


66 Fe 


40 


3.22 


0.51 


68 Fe 


42 


2.30 


0.62 


bu Cr 


36 


1.55 


0.31 


62 Cr 


38 


2.77 


0.66 


64 Cr 


40 


3.41 


0.76 


66 Cr 


42 


2.28 


0.90 



calculation and the one that corresponds to the normal 
filling. The neutron intruder occupations increase with 
/V in both isotopic chains, however the absolute occu- 
pancies are larger in the chromium chain. As mentioned 
already for the N — 40 nuclei, the strong deformation in 
the Cr chain is not only due to the increased population 
of the O59/2 — ld.5/2 doublet, but to the strong proton- 
neutron correlations which tend to be maximal when four 
protons are active in the pf shell as well. 

In Fig. [5] we show also the evolution of the collectivity 
in the nickel chain: in (a) the theoretical energies of the 
first excited 2 + states are compared with experimental 
ones, while in (b) we show the calculated transition rates 
in comparison to the available data. The agreement in 
the calculated energies is very good for all nickels. Con- 
cerning transition rates, the model reproduces well the 
systematics with the minimum in 68 Ni and the rapid in- 
crease of collectivity in 70 Ni. However, the calculated 
value is closer to the lower tip of the error bar. The 
transition rate in 64 Ni is underestimated. In this case 
we know that a better agreement with experiment can 
be obtained in a full pf shell calculation as the neutron 
excitations from O/7/2 orbital are here more important 
than those through the N = 50 gap and, as expected, 
the occupation of the 1^5/2 orbit remains close to zero in 
the calculations. 

Finally, let us note the recent work on 68 Ni [28|, where 
the excited + states have been investigated and a can- 
didate for a proton 2p-2h intruder has been proposed at 
an excitation energy of 2.2 MeV. We have calculated the 
excited + states in our model space and we obtain the 
first excited 0^ state at an energy of 1.2 MeV and the 
second one 0^ at an energy of 2.4 MeV, in good agree- 
ment with experiment. The too low excitation energy 
of the first state can be a result of the delicate mixing 
between the 0^~ ground state and the excited 0^ state. 
These two states are characterized by similar proton oc- 
cupancies with leading Op-Oh (neutron) configuration for 
the Of ground state and 2p-2h (neutron) configurations 
for the 0^. For the 0^ state, located at 2.4 MeV, the 
striking feature is that the dominant proton configura- 
tion has exactly two O/7/2 protons less than the ground 
state, i.e. a pure 2p-2h proton configuration. 

The total quadrupolc sum rules for these states amount 
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FIG. 5: (Color online) Theoretical versus experimental ener- 
gies and transition rates of the nickel isotopes in the vicinity 
of N = 40. 



to 904, 1162 and 2025 e 2 fm 4 for the ground state, the 0^ , 
and the 0^~ states respectively. This shows in particular 
that the 0;j~ state carries moderate deformation, despite 
of its proton intruder nature. 



V. CONCLUSIONS 

The aim of this work was to study the rapid onset of 
collectivity below 68 Ni, suggested by the experimental 
evidence, which is a great challenge for any theoretical 
model. Here we have presented shell model calculations 
in a large valence space including the pf shell for protons 
and the 0/ 5/2 , lp 3 / 2 , lpi/2, 0cf 9 / 2 , and ld 5/2 orbitals for 
neutrons. The effective interaction for this model space 
has been build up using sets of realistic TBME as a start- 
ing point and applying monopole corrections thereafter. 
The progresses in algorithms and computer power have 
made it possible to achieve the largest shell model diag- 
onalizations in this region of nuclei up to date. 

The present calculations have been performed in Ni, 
Cr, Fe, Ti and Ca isotopes around N = 40. We have 
found a satisfactory agreement between the theoretical 
results and the available data, for both excitation ener- 
gies and transition rates. In particular, it has been pos- 
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siblc to describe correctly the rapid onset of collectivity 
below 68 Ni without any collapse of the spherical gaps. 

It has been also shown that the onset of deformation 
develops at N = 40 in the iron chain and already at 
-/V = 38 in chromium isotopes. The maximum deforma- 
tion is predicted for the chromiums which exhibit features 
typical of rotational nuclei. The continuous advances in 
the experimental side, with particular regard to the tran- 
sition probability measurements, will be a stringent test 
for these theoretical predictions. 

The calculated wave functions of the deformed ground 
states were shown to contain large amounts of many 
particle-many hole configurations, with around four neu- 



trons occupying intruder orbitals Ogg/2, 1^5/2 ■ We have 
also related the observed shape change with the under- 
lying evolution of the spherical mean-field, which bears 
many similarities with the one at the N = 20 island of 
inversion. 
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